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Abstract: In this paper we consider the process of diffraction dissociation in deep 
inelastic scattering producing a small mass. This process is analyzed by the calculation of 
qq and qqG production. We show that the small distance contributions (r± oc l/Q) to 
the longitudinal polarised virtual photon dominate the diffractive channels. Formulae for 
the cross section using the gluon density are written within the framework of perturbative 
QCD. We show that the production of small masses by a transverse polarised photon is 
concentrated at moderate values of r_|_ ~ \GeV~ 1 , where the pQCD approach is valid. It 
is shown that only q q pair production contributes to diffraction dissociation with (3 > 0.4 
and the possibility to extract the values of the gluon density from the measurements in this 
kinematic region is discussed. Shadowing corrections are assessed for both longitudinal 
and transverse polarised photons and estimates of the different damping factors are given. 
The relation between diffractive production and the corrections to F 2 is alluded to. The 
evolution of the diffractive structure function is studied and a solution to the diffractive 
evolution equations is proposed. 
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1 Introduction 



Over the last decade, starting with the paper of Bartels and LoewejjTJ, diffractive processes 
in deep inelastic scattering (DIS) have attracted a great deal of attention. The reason 
is that these processes can be calculated within the framework of perturbative QCD 
(pQCD)[Q and their calculated cross sections turn out to be proportional to the square 
of the gluon density xpG(Q 2 ,xp). Q 2 is the virtuality of the photon and xp = ® + , 
where s is the cm. energy squared of the reaction. M is the mass of the diffractively 
produced hadronic system. The recent renewed interest in these channels was triggered by 
Ryskin[|3| and Brodsky et.al.0, who proved that vector meson electro production can be 
calculated in the leading log approximation of pQCD, while the non perturbative effects, 
originating from large distances, can be factored out in terms of the light-cone qq wave 
function of the produced vector meson. 

The goal of this paper is to study single diffraction dissociation (DD) production of 
a small mass hadronic system. We consider a small mass as the hadronized product of 
either a qq pair or a qqG system in the final state of the DD process. We will show that: 

1. This process describes a significant part of the experimental DD cross section. 

2. The small distance contribution (r± oc 1/Q) dominates the longitudinal polarised 
virtual photon which induces DD. This justifies the use of a pQCD approach to calculate 
the relevant cross section. 

3. The production of small masses by a transverse polarised photon are concentrated at 
moderate values of r\ ps 1 GeV~ 2 . This encourages us to believe, that our estimates 
which are obtained in pQCD, could be responsible for a considerable part of the DD cross 
section. 

4. Shadowing (or screening) corrections (SC) are not negligble in the region of low xp 
and should be taken into account. 

The plan of our paper is as follows: In the next section we calculate the process of DD 
by longitudinal and transverse polarized photons into a qq pair as well as through qqG. 
This is done in the ^-representation, which provides a natural frame to assess at what 
distances these processes are viable. It is also a convenient representation to calculate 
SC. Most of the formulae which we obtain for transverse polarised photons have been 
previously derived in the momentum representation, see Refe. || || 0Q. In section 3 we 
extend the formalism suggested by Levin and Ryskin|| and Mueller 0, to the case of a 
penetration of a qq pair with a definite value of r±_ through the proton, as well as through 
the nucleus. We will show that for the case of a diffractive production of small masses by 

*Mark Wiisthoff has done the calculations for a longidudinal polarized photon in his PhD thesis and 
presented the results in the 1995 Durham workshop. 
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longitudinal polarised photons, we obtain an elegant closed expression for the damping 
factor |10| defined by 



p) [wilh SC] 
^[without SC] 



7")2 dt 

" da(r+P^X(M)+p) r.+ h . Qr n l i=0 



where t = —q\. Processes at large transverse distances do not enter the above formula, 
allowing us to calculate the SC within the framework of pQCD. We extend our formalism 
also to DD by a transverse polarized photon. We discuss the large distance contribution 
in this case and show that the SC lead to a decrease of the typical distances that are 
responsible for DD. In section 4 we compare the pQCD predictions with the experimental 
data available from HERA [p|[pl (see also Refs. |I^|H ). We show that pQCD is able 



to account for a considerable part of the observed DD cross section, and discuss which 
distances contribute to the different modes present in the final state of the produced 
system. In all our calculations we use the GRV parametrization [|i~5|j for the gluon density, 
which is the main input ingredient in our work. The GRV parametrization which starts 
from a low value of the photon virtuality, allows us to study the value of typical distances 
in the DD process. In the fifth section we discuss the evolution equations for the DD 
processes. A summary of our results and a discussion are presented in the conclusions. 



2 The diffractive production of qq and qqG systems 
without shadowing corrections 

2.1 Notation 

We list below the notation used in this paper (see Fig.l). 

1. Q 2 is the virtuality of the photon in DIS. 

2. xp = ' , where s denotes the square of the cm. energy and M is the produced 

diffractive mass. Bjorken's scaling variable is given by xb = 

"X R - Q 2 _ £a 

°- r Q 2 + M 2 ~ x P ■ 

4. a 2 = z(\ — z)Q 2 + tuq, where z is the fraction of the photon energy that is carried 
away by a quark with a mass uiq. 

5. k± denotes the transverse momentum of the quark, and f± the transverse distance 
between the quark and the antiquark. 

6. b± is the impact parameter of the reaction and is the conjugate variable to q±, the 
momentum transfer from the incoming proton to the recoiled proton. Note that t = —q\. 
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7. lix denotes the transverse momentum of gluons, which are attached to the quark- 
antiquark pair (see Fig.l). 

8. We use the evolution equations for the parton densities in moment space. For any 
function f(x), we define the moment f(uj) as 

f(u) = I' dx f(x) (2) 
Jo 

Note that the moment variable u is chosen so that the uj = moment measures the 
number of partons and the uj = 1 moment measures their momentum. An alternative 
moment variable N, defined such that N = u + 1, is often found in the literature. 

The x - distribution can be reconstructed by considering the inverse Mellin transform. 
For example, for the gluon density it reads 



xG{Q\x) = — / did x-" g(Q\u) (3) 
2m Jc 

The contour of integration C is taken to the right of all singularities. 

For the solution of the DGLAP equations [IB], as well as the BFKL equation |T7|], one 
has the following form 

g(Q\u) = g{u)e<^ (4) 

where 7(0;) denotes the anomalous dimension which, in the leading log - approximation 
of pQCD, is a function of ^ and can be presented as||18fl 



K ' 71 UJ TV to k uj 5J w 

9. Our amplitude is normalized so that 

and the optical theorem can be written 

a tot = Att lmf(s,0) 
The scattering amplitude in fe^-space is defined by 

a(s,b ± ) = ±-l 4e^ l Si f(s,t = -ql) 



2tt 

In this representation 

a tot = 2 J d 2 b ± Im a(s, b±) 

and 

Pel = J d 2 b ± \a(s,b ± )\ 2 
In what follows we use the notation and normalization of Brodsky et.al.l 



2.2 The cross section for diffractive dissociation 



The cross section for the diffractive process (see Fig.l) 

l*(Q 2 ,x B ) + V - X(M 2 ) + p (6) 

can be written in the form 

da 3a em ^ 2 v - ^ f 1 a f . , 

Xp dx~p~dt ~ 8(2^ \ Zf ^Jo dZ I [ ) 

1/(7* + P - 9(«,*±) + + 9±) + P)| 2 (M 2 + Q 2 ) 5 f M 2 ^ ^ 



2(1-2 



We have assumed in eq.(7) that the small masses are produced by the dissociation of a 
qq system. We will consider later, the corrections resulting from the emission of gluons. 
In eq.(7), Aj denotes the polarisation of the quark, N F the number of quark flavours and 
Zp the fraction of the electron charge carried by the quark with a flavour F. 

Integrating over z and summing over polarizations, we reduce the above equation to 
the form 

da 3o> em o at f i mi d k\ k t 

x riZ-T+ = 7^3 22 Z f N x 



dxpdt (2vr) 3 V J 1 - P M 2 L ~k 



A/ 2 



where z(l — z) = and a = J z(l — z)Q 2 + rrig = Q for massless quarks. The 

region of integration with respect to k±_ is defined k\ < The factor N\ depends on 
the initial polarization of the photon. In the case of longitudital polarization N\= 4. For 
transverse polarization A^a = z 2 + (l — z) 2 . We follow Refs. || || regarding normalization, 
and the definition of the photon wave function. The main task before us is the calculation 
of the amplitude /. 



2.3 The amplitude in the r± representation 

This approach was first formulated in Ref . || and has been carefully developed in Ref . || . 
During it's time of passage through the target, the distance r± between a quark and an 
antiquark can vary by an amount Arj_ oc R E denotes the energy of the pair and 
R the size of the target (see Fig.l). The quark's transverse momentum is k±_ oc — . 
Therefore 

k 

Ar ± oc R-t- < r± (9) 
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holds if 

r 2 ± s > 2m# (10) 

where s = 2mE. 

The above condition can be rewritten in terms of xp 

This means that at small values of xp, the transverse distance between the quark and the 
antiquark is a good degree of freedom |5] P [P| , and the interaction of a virtual photon 



with the target can be written in the form 

/ = / ^ **(r±,k ± ) a(r ± ,ql) W*(Q 2 ,r ± ,k ± ) (12) 

\1/ denotes the wave function of the produced qq pair which is equal to e tr± ' k± . After 
integrating over the azimuthal angle we have Jq( kj_ r±_ ) and Ji( k± r± ) for the longitudinal 
and transverse polarised photon induced reactions respectively. \1/ 7 * denotes the wave 
function of the longitudinally ( L ) or transverse ( T ) polarized photon. These have been 
given in Refs. || |J and are equal to 

^t(Q 2 ,r ±1 z) = Qz(l-z)K (ar ± ) = Q^K (ar ± ) (13) 

and 

n\Q 2 ,r ± ,k ± ) = iaK^arx) ^ (14) 

r± 

&('r±,qj_) denotes the cross section of the qq pair with a transverse separation r±, which 
scatters with transverse momentum qj_. We will discuss it in detail below, considering 
initially the case where q± = 0. 



2.4 cr(r±,q±) at q± = 

The form for a(r±,q 2 _) at q± has been given in Ref.f| (see eq. (H) of this paper), a can 
be expressed in terms of the unintegrated parton density <j), first introduced in the BFKL 
papers |DJ and widely used in Ref . f2Cfl . This function is clearly related to the gluon density 



and can be calculated using the equation 

a s (Q 2 )xG(Q 2 ,x) = [ Q dll a s (l 2 ± ) <j)(l 2 ± ,x) (15) 
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Using the above equation we obtain the result of Ref. ||, which reads 

<^x, ? i) = ^ ^/tfH,*) ( 1 - e"^ ) ^ ^ (16) 



where = ^SSM 

We calculate this integral using eqs.(3) and (4), and integrate over the azimuthal angle. 
Introducing a new variable £ = r_|_/_|_ we can reduce the integral to the form 

, . = 16C7 Ftts r du_ , ,_ 7M r Jt l- mo 



a = ^ si (ri) *' /. rf f kf^t < i7 > 



preforming the integration over d£ (see Ref. ||2 1|| 11.4.18) we have 



a{ruql) . «~ 7M ( ^ r ^-;;:>> , 18) 

In the double leading log approximation of pQCD, 7(0;) <C 1, and we obtain the following 
cross section for Nq = 3 

a{r L ,ql) = ^ vr 2 ri ( xG DGLAP (±x) V (19) 



This result agrees with the expression for the cross section value given in Refs. [Q] and [23 
provided we neglect the factor 4 in the argument of the gluon density. We have checked 
that eq. (2.16) of Ref. j|] also leads to the same result, unlike the value for a quoted in 
the paper (see eq. (2.20) in Ref.H). 



2.5 The amplitude for small mass production at qj_ = induced 
by a longitudinal photon 

In this section we calculate the amplitude for qq production with a mass M. Substituting 
eq. fll3]), in our expression (12) for /, yields 



Making use of the well known properties of the modified Bessel functions we can perform 
the integration. The final general result is 

/ = ^ f *?L g{u) 1 A-yH (21) 



M 2 Nq — 1 Jc27ii yv a 4 y (3 

r (i + 7 H) r (i - 7 (a;)) r (2 - 7 M) r (2 - 7 H) /3 3 ^ ( 2 - 7 M, -1 + 7H, 1, 1 - 

where a = Q 

In the case of the DGLAP approach we can simplify eq. (^1|) as <C 1. Neglecting 

7(0?) with respect to 1, and using eq. (H) we get 

It is interesting to note that the argument of xpG is k\/{l — (3), which means that small 
distances of the order of r\ oc (1-/3) /kj_, contribute to small mass diffractive production, 
especially at (3 — > 1. 



2.6 DD cross section induced by a longitudinal photon at t = 

Subtituting the amplitude in eq. (§) for the cross section, we have for three flavours and 
three colours 

*'s£sl"= W' p3{1 ~ 2pr (23) 

^ ^A; 2 / „, k 2 A 2 1 



xp G{ — , xp ) 



1-/3)' ^; rr^i 



|2 

M 2 



One can conclude from eq. (|23| ) that only small distance processes contribute to the 
cross section. This fact manifests itself in the log integration with respect to k\. Even in 
the double log approximation of pQCD, this integral converges and can be rewritten in 
the form 



1 dx' P 



' 9(x P G( 47^y, x P ) 



V 



dln(l/x' P ) 



The above formula gives the result for the case of a two jet production originating from 
a qq pair. 
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2.7 The DD amplitude at q± = for a transverse photon 



Using the formulae given in the previous subsections, and taking into account eqs.(3) and 
(4) for the gluon structure function, we obtain the following expression for the amplitude 

Making use of the properties of the modified Bessel functions we perform the integration 
over r_|_ and obtain 

8C F a s 2 , f dw A 1 ,a 2 {u) 



> = N^* 4k -]o™ g M*t7 ) (26) 



r(i + 7 (w)) r(i - 7 ( w )) r(2 - 7 ( w )) r(3 - 7 M) _ 3 



r(2- fl ( w ))r(2) 



^^(3-7(^3(42,1-^ 



where a = Q k 



M ■ 



In the DGLAP approach the above equation can be simplified as 7(0;) <C 1, it reduces 
to the form 

2 1QC F a s k\ 3 / fci \ 

; = 7r a^t^ XpG ((T^)'^J (27) 



2.8 DD cross section for a transverse photon into a qq pair 

Substituting the amplitude / in the general formula for the cross section (see eq. (§)), we 
obtain for three flavours 



1 M 2 
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2.9 DD cross section for qqG production 

The emission of one additional gluon is shown in the diagrams of Fig. 2. One can see two 
different classes of diagrams which describe gluon emission. In the first, the emitted gluon 
does not interact with the target as shown in Fig. 2a. The general way to take such gluon 
emission into account, is to use the evolution equations. We shall write down explicitly 
all formulae for the emission of one gluon at the end of this section. The second class 
of diagrams describes the process in which the emitted gluon interacts with the target, 
(see Fig. 2b). These diagrams have to be calculated separately so as to provide the initial 
conditions for the evolution equations [|7|]. 

We start with the calculation of the diagrams of Fig. 2b. We would like to determine 
those contributions where the smallness of as is compensated by a large logarithm. In 
other words, we are looking for contributions of the order of as ln(Q 2 / k 2 ) . Our first 
observation is that such a term does not exist in the case of a longitudinal polarised 
photon induced reaction. Indeed, this fact can easily be seen from the general structure 
of the diagrams of Fig. 2b, namely, they can be written in a general form 

x P -^-\ t=0 (Fig.2b) oc f 1 dz f 1 % I tj± vFV(r ± ,gi = 0,x')[^]* (29) 
dxpdt Jo J XP x' J 2n 

where a is defined by (see eq. fljjj) ) 

«s(4-) / 4 \ 2 

a(r ± ,ql) = vr 2 r 2 xG(-,x) (30) 



with [xG{r 2 ±1 x)] 2 that has been introduced in Ref . ]24| . The explicit expression for (xG) 2 
was derived in Ref.H] and it is equal to 



v2 2 






M r' 4 Jo 



db\[a GG } 2 (31) 



where 



a GG (r>,x>) = ^r' 2 [x' G(^, x') ) (32) 

Substituting eq. flTBI ) for ^\ in eq. (ffTD, after the integration over z and using the prop- 
erties of the modified Bessel functions, we obtain 

roo rf r 2 a S (4-) / A \ 2 

a{r u ql) cx J a vr 2 [xG(^,x)\ (33) 
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Therefore, we have no log contribution from the r± integration. However, if we calculate 
the same integral for the transverse polarised photon we find 

a(r ± ,ql) oc /~ ^ tt 2 ( x)J (34) 

oc r dr' 2 ln(QV 2 ) ( xG(\,x) \ 



Finally, going to the representation, eq. (B3) yields the result of Ref. 



da T ^ A7i 2 a em Z 2 r l dz f — dk 2 



{ (3 2 + (z - (3) 2 } (1 - zf (2z + I) 2 { x P G(k 2 , x P ) 
Taking the integral over z we obtain 



^Z^IhM = E ™ L TT^r- InTTi ( 36 ) 



dx P dt ,l ~ ux ' y Q 2 JQi k* 32tt Ak 

{ ~ ~ — [4 - 31/3 - 63/? 2 + 50/? 3 - Uf3 A ] - \nf3 [1 + 10/3 - 2f3 2 } } (a;p G(x P , fc 2 ) 
6 v 

We now calculate the diagrams of Fig. 2a. For a transverse polarized photon, this has 
been done in Ref. and the result in our notation is 



da T _ An 2 a em Z 2 F f 1 dz dk 2 4a§ , M 2 



* ^ = ? * h 5 4 ^ ln & (37) 

^{/5 2 (z-/3) 2 - //3 2 (l-/3) 2 } (,pG(iVp)) 2 = 
_ 4vr 2 a em Z 2 fi d£; 2 4a| , M 2 

E — ^ — /„ L tt^t ln 



Q 2 7/3 7og A; 4 3AT C 7r 4fc 2 
(3{l - (3) 2 {-2(3 ln/3 + ^L[l + 6/3-9/3 2 - 24/3 2 - 6/3 4 - 4/? 5 ]} (x P G(£; 2 ,xp) 

One can calculate the one gluon emission, for the longitudinal polarised photon, using 
the general formula obtained in Ref. |7[ 

da L ^ A7T 2 a em Z 2 F r l dz 
Xp dx^dt lt=o{Fig - 2a) = E Q2 A, - (38) 
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where $p denotes the usual DGLAP splitting function, and $p is the splitting func- 
tion of the Pomeron into a quark- ant iquark pair. For the transverse polarisation of the 
photon this splitting function has been calculated previously in Refs. f|[[7| and we have 
recalculated it in the r_|_ - representation (see eq. (|28|)) 



l\c 

For the longitudinal polarization we obtain the same factorized answer with one important 
difference, the integration over k 2 is logarithmic (see eq. (|23|)). However, as far as the z - 
dependence is concerned, the formula is the same as eq. ( j38|) with the splitting function 
which we have calculated previously (see eq. (PIQ) and which is equal to 

= |^ 3 (l-2,) 2 (40) 

Finally, the answer for the one gluon emission in the case of a longitudinal polarised 
photon is 

da L . ^ A-K 2 a em Z 2 F f 1 dz dk 2 4ag Q 2 , . 

Xp d^dt^ {Fl9 - 2a) = £ Q 4 h * k * ( } 

P \^{P 3 (z-W ~ z 7 p 3 (l-2{3) 2 }(x P G(k 2 ,xp) 
Ait 2 a em Z 2 F r^r- dk 2 4a§ ,_ Q 2 



V Q 4 -/Qg fc 2 3Ar c vr A; 2 



{-2/3 3 ln/?(l-2/5) 2 + A + i/? + ^ 2 -4/? 3 + 9/? 4 + ^ 5 -^/5 6 -|/5 7 } (x P G(A; 2 ,xp) 



3 Shadowing corrections for the diffractive produc- 
tion of small masses 

The advantage of using the r± representation is not apparent in the calculations performed 
in the previous section, in that we obtain the same results previously obtained by using 
momentum space calculation techniques (see for example Refs. [|] || @ [^4| ) . The above 



calculations are instructive and they provide an explicit example of the relations between 
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different variables, used to describe deep inelastic process in different frames. In this 
section, where we calculate the shadowing corrections for the various processes, the rj_- 
representation simplifies the calculation. We will only calculate the damping factors for 
the penetration of a qq - pair through the target, assuming that all SC to the gluon density 
have been included in the phenomenological gluon density that we use in the calculations. 
The size of the corrections due to damping in the gluon sector has been estimated in 
Ref.|2§. 



3.1 The b± dependance of the amplitude 

To estimate the shadowing corrections, we need to know the profile of the amplitude in 
impact parameter space, which means that we need to know the amplitude not only at 
q±—0, but also at all values of the momentum transfer. The gluon density is weakly 
dependent on q± at small values of q±, both in the double log approximation (see a 



detailed discussion in Ref . pOf ) and in the BFKL approach (see Refs.f27|] p8|1 ). Therefore, 
the leading g^-dependence comes from the form factor of the qq pair with a transverse 
separation of r±, and the form factor of the target (proton). 

As the form factor cannot be treated theoretically in pQCD, we will assume an expo- 
nential parametrization for it, namely 

F p (ql) = (42) 

The slope B can be extracted from the experimental data on hadron-hadron collisions, 
provided we take the Pomeron slope a'=0. Namely, B = B^f(a' = 0), where B e \ is the 
slope in the differential cross section of the proton-proton collision |29|]. It turns out that 
the value of B deduced from the hadronic data is very close to the one obtained from the 
proton electromagnetic form factor. 

The form factor of the qq pair with a transverse separation of r± is equal to 

fw<«i) = n, ( Ckl± = t ] ' ~ ) *£ ~ f 21> " - ) m 

where hi (k'j) denote the momentum of the quark i before and after the collision. Each of 
the wave functions is given by the exponent and a simple sum of different attachments of 
gluon lines to quark lines. Hence, 

2\ ii^A I 1 M 



Fq q (q±) = 1 - e*<— (44) 
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The last factor is absorbed in the expression for the cross section, while the first factor 
gives the q± dependence of the qq form factor, which after integration over the azimuthal 
angle is 

F,M) = (45) 
To proceed with the calculation we require the profile function in b± space 

S(b 2 ± ) = ± J A±e^ F p { q \) F- qq {ql) (46) 
The explicit calculation gives 



,2 



To simplify the calculation we replace the above function by 



S(bl) = -L A (46 



where 



r 2 . 



B ' = B 1 + ~ B 1 + ^> (49 > 



For sufficiently large values of a, we consider a 2 B ^> 1, and neglect the second term in 
our calculation. 



3.2 Penetration of a qq-pair through the target 

To calculate the shadowing correction we follow the procedure suggested in Refs.|HJ. 
Namely, we replace cr(r±, q\ = 0) in eq. ( |T2"D by 

a sc (r ± ) = 2 Jd 2 b ± (l - e -|^x >? i=o)S( 6 i)^ (5Q) 

The above formula is the solution of the s-channel unitarity constraint 

2 Ima(s,b±) = \a(s,b±)\ 2 + G in (s,b±) (51) 

where a denotes the elastic amlitude for qq pair with a transverse separation r±, while 
Gi n is the contribution of the inelastic processes. The inelastic cross section is equal to 

a m = J d 2 b ± G in (sM = \d 2 b L (l - e -^i=°)^i)) (52) 
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Our formulation is based on the physical assumption that the structure of the final state 
is obtained from the uniform parton distribution derived from the QCD evolution equa- 
tion. We neglect the contribution to the inelastic final state of all DD processes with 
large rapidity gap. For example, "fan" diagrams that give an important contribution are 
neglected [ED| as well as DD in the region of small mass, which cannot be presented as the 



decomposition of the qq wave function. 

^From a point of view of Feynman diagrams, eq. fl52"| ) sums all the diagrams of Figs. 
1 and 2 in which the qq -pair rescatters off the target and exchanges many "ladder" 
diagrams, each of which can be represented by the gluon density. This sum has been 
performed by Mueller 0], and we shall perform the calculation for the case of small mass 
DD. 



3.3 Damping factor for diffractive production of small masses 
induced by a longitudinal photon. 



Substituting a sc of eq. (|50D in eq. (|i~2|) and using a in the form of eq. (|19|), enables us to 
estimate the general term of the expansion with respect to a. It has the form 



f = C 2Q J — Ko{ar ^^r- U^r 1 (53) 



n ^ e i>in(i/*p) r(i + 7 (^)) r(i- 7 (wO) 

11 L ^~ 9*W 



2 



c> 2m lJ (r(2- 7 (^))) 

( Jo(k ± r ± ) J d 2 b ± S n (bl) 

We replaced all numerical coefficients by the factor C . 
Integrating over r± and b± we obtain the expression 

r ■■ r 1 ( as4Cp7T V (u) 

} - C 2Q {b'(n% - i) ; (54) 

" r du ie 5:<«wv*p) r(i + 7 ( Wi )) r(i- T ( w< )) 

11 L ^~ 9i{^i 



, JCi 2m ^ l > (r(2- 7 (^))) 2 

2 n / o2 \ n 
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In the double log approximation of pQCD with 7(0^) <C 1, neglecting 7 and taking 
1 — p <C 1, we derive a very simple formula. Taking the integral over Ui yields 



„ (3 2 (-l) n - l n\ ( C F 7i A(3 



,a 



2 



r = cb 'jq ~^r- [ww^j 5 as x ^t Xf) J (55) 

Finally, for / we have 

f = CB S^Q (n ~ 1)! U-l) ^" 5 *pG(j>*p)) ( 56 ) 
The above series can be summed to give the analytic function E±, namely 



1 



where (for N c = 3) 



2 a S 7r(3 a 2 2 a s n (1 - 0) / fcj \ 

Using the above equation, we obtain the following expression for the damping factor 



U k 2 dk 2 E 2 (±)e~« 
D\ = JQ ° g ^ ^ (59) 
IqI k ± dk ± K q 

The behaviour of the damping factor for small and large K q can be found by using the 
well known property of E x . Namely, for n q <C 1, D — > 1, while at K q >> 1 the damping 
factor vanishes as 

M 2 

fn* k \ d k \ In /t 

Di K V 2 (60) 

Jgjj kj_dk^n q 

The behaviour of Df as a function of (3 and Q 2 is given in Fig.3. One can see that the 
value of the damping factor reaches 0.5 in the region of small (3 and Q 2 . This means that 
the SC should be included even for a longitudinal photon, where the smallest distances 
of the order of ^ contributes. 
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3.4 Damping factor for diffractive production of small masses 
induced by a transverse photon 



Repeating the steps of section 3.3 we derive the formula for the damping factor for a 
transverse polarized photon. Substituting a sc of eq. (|50|) in eq. (|12|) we obtain 

f =Ca J —— Kx(ar ± ) - (^T^T J (61) 

n dUif £<*Hi/x P ) r(i + 7 (^)) r(i- 7 (^)) 

\ l k 2m 9i[Ui) (r(2- 7 (^))) 2 

[jj Mk± r ± ) j d% ± S n (bl) 

C stands for all numerical coefficients. 

After integration over r± and 6j_ we have 

f = C ^ U(^-i) J (62) 

n dUi( £<«Wi/*P) r(l + 7 (q;t)) r(l- 7 (^)) ,a\ T M» t ) (FY 

\ [ k 2vr 9i[UJi) (r(2- 7 (^))) 2 V U; 

n n 

r(2+n-^7(^)) r(l+n-X)7(wi)) i*i (2 +n 1 -n + 5>(w*), 2, 1 - /?) 

j i 

Again, in the double log approximation of pQCD with 7(0^) <C 1, neglecting 7 as well 
as considering 1 — j3 <C 1, we take the integral over Ui and get 



^ r w P k ± !±^!M f 9ll ^ r( a2 \\ <G*\ 
f = ° B V n UW-1) ^ ^' } 1 ( } 



Finally we obtain 

The above series can written in terms of the analytic function Ei 

f = CB'^ (l + (1- -)tfx(-) e-) (65) 
2 a 2 V K s K s / 
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where n q is given by eq. (|58|). 



Using the above equation we get 



M 



dk\{\ 



1 - 



E 1 i 



A/ 2 



(66) 



Fig. 4 shows the dependence of D\ on f3 and Q 2 . It is interesting to note that the value 
of D\ is small for small (3. This reflects a well known fact that the SC are large for DD 
producing a system with a large mass pi pO 



3.5 The relationship between SC, F2(Q 2 ,xb) and the DD pro- 
cesses 

The simplest relation between the corrections to F 2 and the DD cross section can be 
derived directly from the AGK cutting rulespl (see ref.flTJ) and it reads 



AF 2 (Q\x B ) 

F 2 (Q 2 ,x B ) a tot 1 ' 

where 

F 2 = F 2 DGLAP - AF 2 (68) 

However, eq. ( |67D is only valid when the DD cross section is small. The notion of what 
is meant by "small" in diffractive production is not unique, since DD is small in two 
cases: i) when the kinematic region between partons is very small, this can be dealt with 
in a perturbative way, or ii) in the case when the interaction is so strong that we have 
scattering off a black disc. As we have shown, the damping factors for DD turn out to be 
rather large. It is therefore necessary to reconsider the simple relation in eq. fl6"7|). 

We generalise eq. ( |5"7| ) by calculating the SC for F 2 due to the penetration of a qq pair 
through the nucleon. As we have mentioned above, we adopt throughout the paper the 
approach that the SC in the gluon sector have been already taken into account in the 



phenomenological set of densities that we use[15|. The expression for the F 2 including SC 



for qq pair has been derived by Mueller and in our notation has the following form 

F 2 {Q\x) = f Q dzj ^± *£(Q 2 ,r x ,z) a sc (r ± ) {Q\ r ± , z))* (69) 
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As was shown in ref. || , whithin the LLA of pQCD which we use throughout this paper, 
we can safely replace ^t^t by -r after integrating over z in eq. (|69|) . Mueller's formula 
finally reads (for Nc = Np = 3) 

F 2 (Q 2 ,x) =Y.*f-J— /T— i (l " exp(-ia(r ± ,gi = 0)5(6i) 

(70) 

Adopting the same procedure of integration as for eq. ( |59f) and eq. fl66f), we derive the 
formula for the elastic damping factor which we define by 

D el = VpNBA\ ( 71 ) 

where F^ BA denotes the correction to F 2 due to two Pomeron exchange. In other words, 
pNBA - g ^ e next order to the Born approximation of eq. ([70|). The formula for D 2 el is 

D 2 = |(l + ^)(£iK) + ln S + C) + 1 - - 2n q \ 

el K 2 \ ) 

T 



where n q is defined by eq. ( |58D at ^ = Q 2 . D 2 el is plotted in Fig. 5 for different values of 
Q 2 . It is interesting to note that D 2 el turns out to be very close to 1 over a wide range of 
xp because D 2 el = 1 — ^ at small values of K q . 

Using D 2 el and D^.{D 2 L ) we can rewrite eq. fl67|) in a more general form 

AF 2 (Q 2 ,x) 



F 2 (Q 2 ,x) 

^ dp da° D (Q 2 ,x,f3) D 2 el + 1 f 1 d dcrE D (Q 2 ,x,0) D 2 el 

&tot(Q 2 ,x) Jo dx P D\ (?tot(Q 2 ,x) Jo dx P D\ 



(73) 



D- 



^From Figs. 3 and 4 one can see that the ratio can be rather large, and reaches a value 
of about 2 at low x in the HERA kinematic region. 

4 Numerics and results 

4.1 What distances are essential in DD processes? 

The first question that we wish to address is, "what are the distances which are important 
in DD processes?" 
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1. In the case of a longitudinal polarised photon, the answer is obvious, just by examining 
eq. (|2"3D . Indeed, in all calculations that lead to eq. (p2|) and eq. (|3|) the typical distances 
are r± oc i = Since the integral over k± in eq. (^J) is logarithmic, we conclude 
that k± ~ M/2, from which we estimate that the dominant distances are r±_ « ^. 
Hence, we can safely utilize the pQCD approach for inclusive production initiated by a 
longitudinal polarised photon. In Fig. 6 we plot our predictions for j3 = 0.8. In spite 
of the l/Q 2 suppression we obtain a large value for the cross section. This suggests an 
alternate way to extract the gluon density, similar to that suggested for vector meson 
production^ |j2"5| . It has several advantages: (i) The cross section is larger than that 
for vector meson production; (ii) The prediction is independent of the form chosen for 
the nonperturbative wave function of the produced vector mesons, which is an inherent 
difficultly when making a definite prediction for the DD cross section for vector meson 
production [|5|] |30) ; (iii) The seperation of the longitudinal from the tranverse polarised 
photon is not difficult. We have not found any contamination fron a transverse polarised 
photon induced reaction, for events with large values of (3 (practically for (3 > 0.7). 

2. In the case of a transverse polarised photon, the value of the distances for which the 
pQCD calculation of the DD cross section is valid, depends crucially on the behaviour of 
the gluon density function at relatively small values of the photon virtuality. Indeed, the 
cross section of eq. (^) has an extra k 2 in the denominator, and the integral on first sight 
appears to be infrared divergent. This is not so, as the gluon density is proportional to 
k 2 //j, 2 at small values of k 2 . This fact is a direct consequence of the gauge invariance of 
QCD. However, there is a danger that at very small values of k 2 in the nonperturbative 
QCD region, we could still have divergencies as /i 2 is about the size of the confiniment 
radius. We know that \i depends on x, and tends to be large at very small x PlpOfl. The 
question arises as to what is the situation in the kinematic region of HERA? In an attempt 
to answer this question we adopt the following strategy: We believe that the available 
parametrizations of the DIS structure functions (GRV|15|, MRS [plR , CTEQ |j33| ) describe 
the behaviour of the gluon density in the region of relatively small virtualities. The 
common feature of all these parametrization is the fact that the behaviour xG(k 2 , x) oc k 2 
starts at virtualities which are sufficiently large in the small x region. For larger x and 
at larger values of k 2 , xG(k 2 ,x) starts to be proportional to k ( xG(x,k 2 ) oc k/jj). For 
example, at Q 2 = 2.5GeV 2 , in the GRV parametrization, xG(x, k 2 ) oc k/fi at x < 3.10 -3 , 
while in the MRS (A) this happens at x < 10~ 2 . Relying on the available sets of parton 
densities (which contain all accumulated experimental information on the subject), we 
conclude that even for a transverse polarized photon, the typical distances in the DD 
processes are not larger than r± « \GeV~ x . Hence, we can apply the pQCD approach 
in this case as well. 
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To illustrate this we plot in Fig. 7 the integrand of eq. (|28|) 



Hk\x P ,P) = ± <- [1-2-^1 , = = (74) 




M 2 



as a function of x P and k 2 at fixed (3. We wish to emphasise that the integrand can 



be viewed as the product of two factors: ^xpG(y^, xp)^ jk\ and the kinematic factor 
(1 — 2 -^j ) , 1 . The first factor has a maximum at k\ ~ 1 GeV 2 for all available 

parametrizations of the gluon density. The origin of this maximum is very simple and 
can easily be traced if one uses the semiclasical parametrisation of the solution to the 
evolution equations in the form 

xpG(k 2 ,x P ) -> (—y^^p) (jfe 2 )-r(*Vp) at Xp o (75) 

Xp 

where both u and 7 are smooth functions of ln(l/xp) and In A; 2 . The common feature of 
the various parametrisations of the gluon density is the fact that 7 > ^ at x P < 10~ 3 . 
Such a behaviour manifests itself in a maximum of the first factor at k 2 m lGeV 2 . It 
should be stressed that the argument of the gluon density in eq. fl28|), namely k\j (1 — (3), 
leads to a substantial increase of the typical transverse momentum in the integral. This 
amounts to a decrease of the typical value of the distances important for diffractive q q 
production. The kinematic factor also tends to increase at large values of k 2 . 

Finally, we would like to repeat our statement that the DD process initiated by either 
longitudinal or transverse polarised photons occur at small distances. This fact is essential 
to justify the use of pQCD for the calculation of DD. 



4.2 The xp dependence, factorization and the "Pomeron struc- 
ture function" 

The xp dependence of the calculated F^ D ^ is plotted in Fig.8 together with the HI 
data [TJIj, where 

x P F2 D{3 \Q 2 ,xp,f3) = —Q f dt x P da = ® x P f a , \ t =o (76) 



An 2 a em J dxpdt ABn 2 a ern dxpdt 

and B denotes the slope of the DD cross section. We take B = 4.5 GeV~ 2 in accord with 
the the prelimenary experimental data from HERA (see Ref. [35|). The value of B has 
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not yet been measured with good accuracy. The above value of B coincide with the slope 
that has been extracted from the high energy phenomenology of " soft" processes (see, for 
example, Ref.ES]). 



We wish to stress that in all our comparisons with the experimental data, including 
Fig. 8, we include the SC, which have been calculated in the previous section. This means 
that the formulae for the cross sections, obtained in section 2, are multiplied by the 
damping factor D\ or D\, calculated in section 3. 

We now concentrate on the xp dependence of F® D ^ . Based on the Ingelman-Schlein 
hypothesis of the Pomeron structure function [34], F® D can be written in a factorised 
form 

F 2 DD{3) = fixp) Ff(Q 2 ,/3) (77) 

where F 2 P denotes the Pomeron structure function and f(xp) the Pomeron flux factor. If 
the Pomeron is a simple Regge pole then f(xp) oc (^) 1+n - n is intimately related to 
the intercept of the Pomeron trajectory: ap = a p (0) + a' P \t\, where op(0) = 1 + e. 
Thus: n = 2ap(0) - 2 = 2e. 

Our approach is orthogonal to the IS one. We do not expect eq. ( ff7|) to scale, and if 
we use the parametrisation suggested in eq. ([77]), and fit the xp dependence assuming the 
form 

F DD(Z) = ( J_ r F P {Q 2 ji3) (7g) 
Xp 

we expect n to depend on both (3 and Q 2 . 

We present the calculated values of n(Q 2 ,/3) in Fig. 9a. We would like to point out 
that the value of n is much bigger (n rs 0.55) than the experimentally fitted value of 



(n ~ 0.2) (see Ref.[14]]). The difference is due to the steep behaviour of the gluon density 
function in the GRV parametrisation. Notice, however, that the new ZEUS data||13|| 
gives a value of n w 0.4 which is much closer to our result. Direct comparison with 
the experimental data given in Fig. 9a, shows that our calculation is able to describe the 
experimental data at xp ~ 10~ 3 but it substantially overshoots the experimental data for 
smaller values of xp. We think this behaviour is an artifact of the GRV parametrisation 
in the leading order of pQCD. A simple ad hoc procedure to improve our results would 
be to take GRV in the next to leading order. In our opinion such a procedure is not 
self consistent unless the entire calculation is done in a higher order. One can see from 
Fig. 9b, where the ratio of the leading order GRV gluon density to the next to leading order 
distribution (x P G(Q 2 , x P ) LO /x P G(Q 2 , x P ) NLO ) is plotted that the use of x P G(Q 2 , x P ) NLO 
instead of xpG(Q 2 , xp) LO suppresses the value of the diffractive cross section at small Xp. 
It is interesting to note that the effective power n reduces then to a value n w 0.3 
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which is closer to the experimental value suggested by the new ZEUS data[[L3|]. The 
second source of suppression at low xp comes from the energy behaviour of the slope 
B, where B = B + 2a' P \n(l/xp). Note that for a Regge soft Pomeron we have 
otpif) = 1 + e + a'pt. 

Regardless of the above, the fact that pQCD calculations predict a larger cross section 
than the experimental one, lends support to our statement in the previous subsection 
that at small distances, where we can trust the pQCD approach, the calculation of the 
diffractive production of small mass, appears to work. How the xp dependence is affected 
by the different parametrizations of the structure function is still an open question. We 
plan on clarifying this in a future publication. 

4.3 The diffractive structure function F® D ^ 
In Fig. 10 we display our calculated values for F 2 DD ^ defined as 

F 2 DD(2) = [ (XPUaX dx P F 2 DD(3) (79) 

where the values of (xp) m j n and (xp) max are taken to have the same values as in the 
relevant experiments. In Fig. 11 we compare the calculated F 2 DD ^ values with the relevant 



HI data|14|. /,From these figures we are able to conclude: 

1. The pQCD calculation provides a fair description of the experimental data. 

2. The contribution of the longitudinal polarised photon is important at large (3 > 0.7. 

3. The process with an emission of an extra gluon gives rise to a sizeble contribution for 
(3 < 0.4 and should be taken into account in a more consistent manner by solving the 
evolution equations for DD processes (see the next section). 

4. For (3 > 0.4, the emission of an extra gluon provides only a small contribution, so one 
can attempt to extract the gluon density by measuring F 2 DD ^ at large (3. 



4.4 The transverse momentum spectra 

Our formulae ,see eqs. (23), (28), (36), (37) and (41), can be used for a more detailed 
analysis of the DD events. In particular, we can describe the transverse momentum 
spectra of the parton jets in DD. Basically, these spectra are given by the integrand of 
our formulae, as the kj_ which appears is just the transverse momentum of the produced 
jet. 
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K ~ ^dW) ^ 



In Fig. 12, the ratio 

F 2 DD{2 \kl > kl) 

is presented with different values of h\. The physical meaning of this ratio is that it 
indicates the fraction of the all DD events possessing transverse momenta larger than 
ko- One can see from Fig. 12 that we expect a fairly large fraction of the events with big 



transverse momenta, this is in agreeement with the new HI data|36 
In Fig. 13 we plot the calculated value 

XpF DD(i) _ Q 2 x P d 3 a | 
An 2 B a em dxpdtdkj_ 

where k± denotes the transverse momentum of the jet. The comparison with the ex- 
perimental data is also quite good, at least, the k± distribution reproduces the main 
qualitative feature of the experimental data |36| , namely, smoother behaviour than at 



small values of k\, and tj- behaviour for larger k±. On the other hand, we predict a sharp 
drop at large k±, which has not been seen in the data as yet. Perhaps, the emission of 
two gluons starts to be important at such large values of transverse momenta. 



4.5 Matching with diffractive vector meson electroproduction 

Vector meson diffractive electroproduction corresponds to the kinematic region of (3 — > 1, 
namely 

a Q2 - i M2y rsn 

where My is the mass of the produced vector meson. In this kinematic region we can 
rewrite eq. fl23D and eq. fl28D in the form 



da L 

dAPdt = (82 > 
i^a em al 2 Ml ,4- dk\ 1 I „(k 2 Q* \V 

V w 

1 da T , 
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^ 2 a em a 2 v2 M v f4 dkj 1 fcl / r (k 2 Q 2 W 

These formulae give the cross sections for production of all hadrons with mass My, and 



reproduce the main features of the exclusive vector meson production||[[|]p5]]. It is inter- 



esting to calculate the cross sections in the region of My = m where we do not expect 
any other mesons, besides the p-meson to be produced. In the double log approximation 
of pQCD, we can put the argument of the gluon density equal to Q 2 . Evaluating the 
remaining integrals we obtain 



° min 



a 



T Mv 



v 



The value we take for Q 2 min depends on how strong our belief in hadron-parton duality is. 
If we believe that the hadron-parton duality can be used at large distances, we can choose 
the low limit of integration from the condition that the argument of the gluon density in 
eq. ([82] ) and in eq. ([83]) is bigger than the intitial virtuality in the GRV parametrization, 
i.e. Q 2 = Q 2 min Q 2 /M v or Ql m = Q 2 M V /Q\ For such a value of Q 2 min we have a L /a T = 
ln^ 

T7T m -753-. 



Experience gained from e + e— annihilation indicates that Q min could be very small, 
about m^. [p7| . In this case, we still have the Q 2 rise of the ratio, but the coefficient is very 
small, a L jo T = 0.125Q 2 . In Fig. 14 we plot o L and a T , using eq. (|82]) and eq. fl8~3|) with 
a cutoff Q 2 min = ml. 

Concluding this section, we would like to reiterate our claim that the simple pQCD 
calculation provides a reliable basis for the discussion of the origin of DD in DIS. It also 
suggests a new way to determine the gluon density. This new measurement can be done 
in two unambigous ways: 

1. the measurement of F^ D ^ with j3 > 0.7, allows us to extract the longitudinal polarised 
photon structure function which is only sensitive to small distances of the order r± m 

i/Q- 

2. the measurement of F® D( - 2 ^ with k± > k > 1 GeV 2 . In this case only distances 
r±_ < 1/ko contribute to our formulae, and we have reliable pQCD predictions for the 
fraction of the DD events. 



5 Evolution equations for DD 

In this section we discuss the evolution equations for the DD structure functions that 
have been introduced in the previous section. The evolution equations for the transverse 
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polarised photon have been proposed and derived in Ref.[[7|, we will comment on them 
later. We first discuss the evolution equations for the longitudinal polarised photon, which 
have not been formulated previously. These provide an interesting theoretical insight into 
the structure of the evolution of exclusive processes in QCD. 



5.1 Evolution equation for DD of a qq-pa\r production by a lon- 
gitudinal photon 



We start with the evolution equation for the simple case of qq production (see eq. (p3|)). 
It is useful to introduce the DD structure function^] absorbing all extra powers of Q 2 in 
the definition 



ol% ^ f^r dk 2 , ( „/ k 2 , x x ~ 



Taking the derivative with respect to InQ 2 , we obtain the evolution equation 

dF° D (Q 2 ,x P ,P) a 2 2 2 / (Q 2 \\ 2 

d^Q 2 - (1-/3) (1-2/3) [vpG^xp)) (86) 

We first investigate eq. fl5E| ) by considering the contribution of the diffractive processes 
to the deep inelastic structure function (i.e. to the total cross section). We integrate 
eq. fl86D as well as eq. (^) over all values ofl > xp > xp. The solution of the 
evolution equation for the integrated F£ D (Q 2 ,f3) {Fgg t = £ g dx P Fg D (Q 2 , x P , p)) 
can be obtained by transforming to the moment representation with respect to \o.{1/xb) 
(see eq. (0)). This has the form 



d\nQ 2 



o^fSntiQ 2 ^) (87) 



where we denote by f G2 (Q 2 ,u) the moment of the right hand side of eq. ( |86| ) inte- 
grated over xp. Our calculation is done in the double log approximation where we know 



that rap? 

j:G 2 trtl \ ,G 2 /^2 \ T2(^)lng2 
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and 

4N c a s 1 A^ c a s 1 , . 

72(0;) = > 7i(^J = 7(<^J = (89) 

H UJ H UJ 



The solution of eq. (p7|) is 

/L;int(Q M = «s A — e Q o (90) 

72 J 

For the deep inelastic structure function, the process of DD initiated by a longitudinal 
polarized photon contributes to the twist four term in the Wilson Operetor Product 
Expansion (note the additional factor of in eq. (0)). We claim that the anomalous 
dimension of the twist four operator is equal to 72. It should be stressed that the diffractive 
cross section does not describe the complete twist four contribution to the structure 
function. The additional gluon emission which increases the value of the total anomalous 



dimension 38 does not contribute to the diffraction into small masses. 



5.2 Evolution of the DD structure function (greneral consider- 
ations) 



In the general case we can rewrite eq. ( pop in the form 
pDD _ «§ f Q2 ^ dkl p 3 s (3 Q 2 F ( ( kj \\ 2 

where £ denotes the singlet structure function of the quark-antiquark pair (q(x, Q 2 ) + 
q(x,Q 2 )) and the Pomeron splitting function <3>£ has been defined above (see eq. (f40"D). 
It should be stressed that only S| enters the evolution equation in the leading order of 
pQCD. 

Differentiating the above equation we obtain the evolution equations for the monents 
in the form 

d\nQ 2 ~ (92) 

where P are the kernels of the DGLAP evolution equations. We obtain the ordinary 
DGLAP evolution equations with an inhomogenious term. This equation has the solution 

fESMi") = ( 93 ) 
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a S rG 2 /^2 w 1 1 t 72 M In §7 

JV C 7 a (w) - 7+(^) gf 2 (w) - 7 (w) 



7+(w) - 7 (w) 7 + M - 7 H 

where the functions fEiLtiQl^ 00 ) an d f G2 (Qo, denote the intial conditions for the evo- 
lution equations, and describes the large distance physics. It is obvious that the second 
term of eq. (p3|) is the solution of the DGLAP equations for fixed as, which we assumed 



in order to simplify the solution. One can find this form of the solution in Refs. |39] |40 



and in the references contained therein. 7 + ' denotes the eigenvalues of the anomalous 
dimensions matrix and in the leading order of pQCD they are equal to 

1 



where 



7* = ^{P§ + P§ ± \/(P§ - Pa) 2 + IPSPn (94) 

r,r «s r „ Ar r 11 1 2 1 1 2TpNp 

po = -^-{2N C [- + - - — - + — — - — — _^(a, + 1+^(1 " — ^1 

P* = ^l {--— + —} (95) 
pS = 2T R N F a s J2 2_ _J_ 

It should be stressed that in pQCD only 7+ has a singularity in the second term at small 
lu, and 7 + <C #2(^0 at w — >• 0. This means that only the first term survives at small x. 

We do not claim that this result is valid at small virtualities, where pQCD arguments 
are not applicable. The value of the DD structure function at a small vituality of the 
photon (Q 2 = Qq) which can be measured experimentally, defines only the specific sum 
of singularities, namely 

f££?(<? ,*B) = (96) 

J- / du (^f'iQi^i 1 + 1 1 + fr s (Qi^)) 

2m Jc \N C 7 2 (w) - 7+ (w) 72M - 7 M / 

The fact that these two terms have different anomalous dimensions, suggests that one 
might be able to use their xp dependence in order to distinguish between them. 

The second term is proportional to ( xpG(Ql, xp)) 2 for the unintegrated F DD , but 
there is no theoretical argument why such a behaviour should be valid for a small value 
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of Qq. Phenomenologically, we associate these two contributions with the "soft" and 
"hard" Pomeron. This terminology reflects our hope that their xp behaviour is different. 
Assuming that their behaviour can be approximately parameterized as x P e , we expect for 
the "soft" (second) term e so f t ~ 0.08 - 0.1, while for the "hard" one (first term in eq. (|93|)) 
(-hard ~ 0.2-0.5. We wish to stress that the problem of separating these two contributions 
is beyond the scope of pQCD, and has to be tackled using a nonperturbative approach. 
As a guide to phenomenological applications we would like to mention that the Ingelman 
and Schlein approach ]33] based on the Pomeron structure function, means that only the 
"soft" contribution survives, while the "hard" Pomeron approach deals only with the first 
tern in eq. (|93|). Returning to F DD (see eq. fl9T|)), we can write the evolution equation at 
fixed xp 

9 -l 1 iz Pl (z) F^, X J-) + (97) 



(9* 



d InQ 2 Jp 4W A ^ ' c ' z- 

f p dzPS{z)F^ G {Q\x P ^) + ^ /3 2 (1 - 2/3) 2 (x p g(^,x p ^J 
In moment space with respect to ln(l//3), the equation reduces to the form 

dfE DiS (Q 2 ,x P ,u) 
dlnQ 2 

P S s{uj)f^ S {Q\x Pl u) + P£(uj)f? D -> G (Q 2 ,xp,uj) + T(w) (x p G(Q 2 ,xp) 
where 

TE;M = - + -^T - Uo \ (99) 

u + 2 uo + 6 uo + 4 o 

The general solution of the above equation is 

fElMM = (ioo) 

JMc 7 — 7 {uo) 7 + {oj) — 7 (cj) 

+ TT-^ + — ^ T7-7 / . T^ T H \XpG(Q ,X P ) ) 

72 (w) - 7 + M 72 (w) - 7 M •'Q 2 , Q v 7 

— / DD _i_ t DD 

— J SOFT i J HARD 

The first term FgQ FT is associated with soft diffraction, and it requires the phenomeno- 
logical input f DD (uo,xp,Ql). We choose for this input 

f DD (Ql,x P ,uo) = T(uo) (xpG(Q' 2 ,x P )) 2 (101) 
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We still need to fix the value of Qo in our initial condition. There is no reason to assume the 
same value of Qq in DD and DIS. We choose, for DD, a value of Qq which is approximately 
1—2 GeV 2 . In doing so we can use pQCD to evaluate the two terms appearing in eq. ( [LOUp , 
and this leads to eq. ( |101| ) for the first term. 

In the /^-representation the intial condition is simply 

FE D (Ql,x P ,(3) = |M 2 (1 - 2/?) 2 (ipg(|x p )) 2 (102) 
This is plotted in Fig. 15 for different values of xp. 



5.3 Evolution equation for transverse polarized photon initiated 
DD 

The evolution equations for a transverse polarised photon have been given in Ref.[^]. 
Introducing the densities qp(Q 2 ,x) and Gp(Q 2 ,x) which denote the quark and gluon 
densities in the Pomeron, we can rewrite the evolution equations in the form (see ref. 
for details): 

dS P (Q 2 ,x) _ 

d\nQ 2 ~ [W6) 

J 1 ^P#(z)S P (g 2 ,^) + P£(z)G P (Q 2 ^) + ^^p-AP) (x P G(Qlxp 

dG P (Q 2 ,x) fi dz 



d In Q 2 Jp z 

{P°{z)^p((l-{3)Q 2 ^ + Pg(z)G P ((l-P)Q\-)} 

+ ^2®%t(P) (x P G(Qlxp) 

where £ = q(Q 2 ,x) + q(Q 2 ,x). 

On inspection we note that the structure of the equations is the same as for the 



longitudinal polarised photon, eq. ( |103| ), and is a normal DGLAP evolution equation 
with an inhomogenious term. P are the kernels of the DGLAP equations while $ denotes 
the Pomeron splitting functions 



30 



G _ a s N 2 1 2 2 

$P;T " 2(^-1) * (i Zj + 

The set of eq. ( |103| ) can be solved in the moment representation using the standard tech- 
niques (see Refs. ||39|||[40|| and references therein). The main difference in the comparison 
with the lomgitudinal photon, is the fact that the inhomogenious term has an extra 1/Q 2 
suppression, provided we choose the appropriate initial condition for eq. (|103|) . The con- 
tribution of this term to the evolution equation is consequently small. To understand the 
structure of the result we write the solution of eq. ( |103| ) in the region of small (3 (oj — > 
), where the emission of extra gluons and quarks could be important (as have been seen 
in section 3). In the region of small oj we can neglect 7~ (see for example Ref. |[40|| ) and 
the solution has a simple form 

Z p (Q 2 ,oj) = (105) 

{ H P {Qlu) + G P (Qlu) Yfai") }e7+ (.)MQVQg) 

^(oj) — 7 (oj) 7 + (w) — 7 (oj) 



+ /, ~7^i ="ri 171 {xpG{Q ,xp)\ 

where g-iiijj) — > 4Ncas - 1 at w -> 0, and $p(cj) is the oj image of the Pomeron splitting 
functions of eq. ( |104p . They are of the form 



No uj + 2 oj + 3 oj + 4 

P;T ~ 2(N 2 - 1) ^ + cj + 1 ~ + 2 ~ oj + 3 ~ oj + A + oj + 5 * 

To fix the two initial moments Ep(Qq,c<j) and Gp(Qq,oj) which enter the solution of the 
phenomenological term fsoFT a ^ l° w ft> we need to choose a phenomenological input. We 
use Regge phenomenology to fix this termf4l| rewriting it in the form 



Xp 

where F dd (Q 2 ,oj) is the oj -image of G PRP /3 A «(°>(1 - (3) n . Ap(0) = a R (0) - 1 where 
or(0) is the intercept of the secondary Regge trajectory (ctp(0) ~ 0.5). The power n 
is not defined in Regge phenomenology and we only know that n > 1. Gppp denotes 
the triple Pomeron- Regge- Pomeron vertex |4T|. The second term is the hard diffractive 
contribution, which is calculated theoretically. 

There is an alternate way to implement the initial condition, i.e. to use the results 
of our calculations in the low mass region. We showed in our calculations for small mass 
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production, that the typical virtuality is sufficiently large. Therefore, we can choose the 
value of the initial virtuality Qq for the DD process to be about 1 — 2 GeV 2 and calculate 
fsoFT m pQCD. In this case 



Zp(QIoo) = -^$% T (u) (x p G{Q'\xp)) 2 (108) 

VO 

and 

G P (Ql,cu) = <f>$(u) (x P G{Q'\x P )) 2 (109) 
This input corresponds to the initial distributions in the j3 - representation 

Zp(Ql,x P ,P) = ^$| ;T (/3) (xpG(Qlxp)) 2 (110) 

VO 

G P (Ql,x P ,P) = (x P G(Ql,xp)) 2 

VO 

These distributions are plotted in Fig. 16 for different values of xp at two Qq values 
(Ql = 1 GeV 2 and Q 2 = 2 GeV 2 ). 

We wish to emphasis two points. First, we propose to seperate between the soft and 



hard diffraction processes in a way which conforms with the factorization theorem p2 
for the integrated diffractive structure function. Although, there is no proof that the 
factorization theorem holds for the case of diffractive contribution, in DIS|43|], all known 



contributions which violate the factorization theorem turn out to be small J38 1 f44| f45 



The second point concerns the value of the scale Qq. There is no reason to take 
the same value of Qq, both in the evolution equation and in the DIS structure function. 
Depending on the choice of the scale Qq, we can start the evolution with very small 
values of Qq and use Regge phenomenology as an initial condition for the DD structure 
function. Or, we can start with a sufficiently large value of Qq in the DD processes, and 
reconstruct the initial distributions using pQCD. We followed the second alternative, and 
our calculations for small masses supports this strategy. 



6 Conclusions 

DD processes in DIS provide a new window to collective phenomena in the parton cascades 
since they originate from parton-parton interactions and vanish if there is no interaction 
between the partons. However, DD processes possess an inherent awkwardness in that it 
is difficult to predict what is small or large in DD. Indeed, in two extreme limits for small 
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and for large ("black disc") parton-parton interactions, the cross section of the DD turns 
out to be small. In this paper an attempt is made to discuss the DD process within a 
pQCD framework. 

Our main results are: 

1. The DD cross section originates from small distances (r_i_ < 0.2 fm) even for a transverse 
polarised photon. This fact justifies our use of pQCD for DD processes. 

2. For (3 > 0.4 only production of qq pairs contributes to the cross section of DD. This 
suggests a method to extract the value of the gluon density from the DD measurements 
in this kinematic region. 

3. For (3 > 0.7 the longitudinal polarised photon provides the dominant contibution 
to DD processes. DD of a longitudinal photon is one of the most promising processes 
for pQCD calculations, since only small distances (r± oc 1/Q) contibute to this process. 
Therefore, the measurement of DD at such large values of (3 allows one to extract the value 
of the gluon density. It has all the advantages of vector meson production in DDj||[25 



but without any of the uncertainties due to our poor knowledge of the hadronic wave 



function [[ftj [[46 



4. The general approach to the evolution equations for DD has been formulated and 
solutions to the evolution equations have been found. The main difference between the 
evolution equations for DD and the DGLAP evolution equations for the DIS structure 
functions, is the appearance of the inhomogenious term in the DD equations. The form of 
this term and its influence on the solution has been discussed. The inhomogenious term 
in the DD evolution equations (eq. (103)) changes the behaviourof ®^q 2 at large values 
of (3 adding a positive contribution. This contribution may change our conclusion that 
the Pomeron gluon density peaks at /3 — > 1, which was based on the DGLAP evolution 
without the inhomogenious term. 

5. The formulae for SC have been obtained and the damping factors have been calculated. 
The main result is that SC are important in the case of DD, and the value of the cross 
section for DD processes allows us to estimate the corrections to the DIS structure function 
which turns out to be fairly large, -^P- ~ 2 ^— . 

We consider our paper as a first step in the understanding of DD as a new source of 
information about collective phenomena in DIS in the region of small x. We hope this 
paper will give some impetus to treat DD as a "hard" process, within the solid theoretical 
framework of pQCD. 
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Figure Captions. 



Fig.l: Diffractive production of a quark- ant iquark pair. 

Fig. 2: Extra gluon emission in diffractive production. 

Fig. 3: The damping factor (D\) for a longitudinal polarised photon. 

Fig. 4: The damping factor (D^) for a transverse polarised photon. 

Fig. 5: The damping factor (D 2 el ) for F 2 (Q 2 ,x). 

Fig. 6: DD by longitudinal and transverse polarised photons 

at = 0.8. 

Fig. 7: The integrand of eq. (28). 

Fig. 8: x P F 2 DD ^ versus x P , calculations and HERA data. 

Fig. 9: Effective power n(Q 2 ,{3) and the ratio xpG(Q 2 ,x p P ) N ° m ^ e GRV paramerization. 

Fig. 10: Different contributions to F® D ^ 2 \ 

Fig. 11: Comparison of _p 2 DD ( 2 ) with the HERA experimental data. 

Fig. 12: The ratio R (eq. (80)) for different values of k . 

Fig. 13: The transverse momentum distribution for the DD jet. 

Fig. 14: The cross section for p production (in arbitrary units) for eqs. (82) and (83). 

Fig. 15: The initial parton distribution for the DD evolution equations with longitudinal 
polarised photons. 

Fig. 16: The initial parton distributions for the DD evolution equations with transverse 
polarised photons. 
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Figure 1: The diffractive production of quark - antiquark pair. 



Fig. 2a 




Fig. 2b 

Figure 2: The extra gluon emission in the diffractive production. 




Figure 3: The damping factor for longitudinally polarised photon. 




Figure 4: The damping factor for the transversely polarised photon. 




Figure 5: The damping factor (D 2 el ) for F 2 (x,Q 2 ). 
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6: The DD by longitudinally and transversely polarised photon at /3 



= 0.8. 




Figure 7: The integrand of Eq.(28). 




Figure 8: xp F 2 versus xp. 




Figure 10: Different contributions to F® D ^ 2 \ 
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Figure 11: Comparison of F 2 with the HERA experimental data. 




Figure 12: Ratio R (Eq.(75)) for different values of k . 




Figure 13: The transverse momentum distribution for the jet in the DD. 




Figure 14: The cross section for p production (in arbitrary units) for Eqs.(77) 
and (78). 
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Figure 15: The initial quark distribution for the evolution equation for lon- 
gitudinally polarised photon. 




Figure 16: The initial quark distribution for the evolution equation for trans- 
verselly polarised photon. 



